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We study the topological structure of matter-light excitations, so called polaritons, in a quantum 
spin Hall insulator coupled to photonic cavity modes. We identify a topological invariant in the 
presence of time reversal (TR) symmetry, and demonstrate the existence of a TR-invariant topo¬ 
logical phase. We find protected helical edge states with energies below the lower polariton branch 
and characteristic uncoupled excitonic states, both detectable by optical techniques. Applying a 
Zeeman field allows us to relate the topological index to the double coverage of the Bloch sphere by 
the polaritonic pseudospin. 


Recently, topologically nontrivial states of matter with 
protected edge or surface states have attracted much at¬ 
tention [1 -3]. While the first realizations were found in 
electronic systems [4-7], topologically nontrivial phases 
of periodically driven systems [8-11] and photons [12-17] 
have been discovered within the last years. A time re¬ 
versal (TR) invariant topological phase can exists when 
a band inversion occurs, as a function of momentum, be¬ 
tween orbital states with different parity, and when these 
orbital states are coupled by spin-orbit interaction. 

In electronic systems, the time reversal operator 
squares to minus one, T 2 = —1, implying the existence of 
degenerate Kramers pairs, such that the crossing of topo¬ 
logical edge states is protected. In contrast, for bosonic 
systems with T 2 = +1, in general there is no TR invari¬ 
ant topologically nontrivial phase [18-20]. The photonic 
topological insulators [12-17] either break TR or have a 
built-in degeneracy, which protects edge states in a way 
similar to the Kramer’s degeneracy in fermionic systems. 

We consider strongly coupled light-matter systems 
in two-dimensions, so-called polaritons, in which the 
bosonic polariton can inherit its topological properties 
from the electronic part. Building on such an example, 
realized by quantum spin Hall (QSH) electrons [4, 6] cou¬ 
pled to cavity photons, we develop a framework which 
allows to characterize topological states of polaritons (cf. 
Fig. 1). In contrast to recent proposals of topological po¬ 
laritons [21-23] where TR broken systems were discussed, 
we focus on TR invariant topological polaritons. Here, 
we go considerably beyond previous works and i) define a 
topological invariant for TR symmetric bosonic systems, 
ii) explain that contrary to [18—20] a topologically non¬ 
trivial phase is possible due to a vortex-like singularity in 
the exciton-photon coupling, and iii) describe how TR- 
invariant topological polaritons can be detected experi¬ 
mentally by looking for dark excitonic states and edge 
states, and by studying the polarization in the presence 
of an external Zeeman field. 

Main results - A polariton consists of an exciton cou¬ 
pled to a cavity photon, such that lower (LP) and upper 
polariton (UP) dispersion branches are formed [24, 25]. 



Figure 1. (color online) Defining properties of a topologi¬ 
cally nontrivial polariton: a) the LP and UP branches (orange 
lines) are split by spin-orbit coupling, except at TR invariant 
momenta (T,Q) (thick black dots). For topological polari¬ 
tons purely excitonic states (gray dot) emerge along a line 
q = qo, and helical edge states (thick purple line) are present 
below the LP. b) with a TR breaking Zeeman field Bz ^ 0, 
the polarization vector n of each LP-dispersion branch can be 
tracked over the whole Brillouin zone. The polariton is topo¬ 
logically nontrivial if n covers the entire Bloch sphere (inset). 
Here, the north (south) pole (blue (red) color) stands for right 
(left) circular polarization, while the x-y plane represents lin¬ 
early polarized light. 


The two polarization directions of the photonic compo¬ 
nent can be identified with a pseudospin [26], which can 
be described by an effective Hamiltonian, and whose di¬ 
rection can be observed by detecting photonic emission. 
For the LP branch it takes the form 

Hlp(q) = eLp(g) + a-h(q) , (1) 

where h(q) is a momentum dependent effective magnetic 
field. With clp (q) we denote the lower polariton disper¬ 
sion in the absence of pseudospin coupling, and u is the 
vector of Pauli matrices. TR invariance TFf(q)T _1 = 
H(—q) with TR operator T = — <J X K. requires that the 
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x- and y-component of the effective magnetic field h are 
even functions and the 2 -component is an odd function 
of momentum. The pseudospin polarization of the eigen¬ 
state Jxi 2 ) is given by n li2 = <Xi,2 |ct|xi,2> = + h with 
h = h/\h\. If n points north (south) the emitted light 
is right (left) polarized, while the x-y plane represents 
linearly polarized light, see inset Fig. lb. In the presence 
of an additional parity symmetry, Hi,p(—q) = Hpp^q), 
one finds h z = 0, and the emitted light is always linearly 
polarized. Physical mechanisms which lead to a non-zero 
effective magnetic field are, for example, a longitudinal- 
transverse splitting of the electromagnetic field [26] or a 
spin-orbit coupling of the electronic building blocks as 
discussed below. In both cases, it turns out that 711,2 
winds twice around the z-axis if q encircles the T-point. 
Due to continuity of h(q), we find that h(T) = 0, implying 
that the LP (UP) dispersion is degenerate at the T-point, 
see Fig. la; in our QSH-polariton model, h(T) = 0 at all 
other TR invariant momenta as well. 

The Chern number [27, 28] for the eigenstate |xi, 2 ) 
counts how many times ni >2 wraps around the unit sphere 
if q covers the Brillouin zone, C 1,2 ~ + J- h■ (d qx h x d qy h). 
Clearly, in the presence of parity symmetry, Ci i2 = 0 as 
h z = 0. However, even without parity, Cp 2 = 0: The 
scalar triple product contains each component h z with 
i = x,y, z exactly once, and thus is an odd function of 
momentum, because h x>y is even and h. z is odd. Since 
the integral over the Brillouin zone is invariant under 
momentum inversion, the Chern number has to vanish, 
in agreement with the general classification [18-20]. 

For polaritons in a QSH cavity the LP eigenstates are 
superpositions of excitonic |6i, 2 ) and photonic wave func¬ 
tions |a lj2 ): |$i, 2 ) = /3 i,2 1 ^ 1 , 2 ) +au, 2 |ai, 2 ) with real co¬ 
efficients /3i,2, cki ; 2, see below for details. An explicit 
calculation shows Tj$ lj2 (g)) = = F| < f > i, 2 (—?))■ By con¬ 
struction, the effective model is the projection onto the 
photonic sector, such that |,\i, 2 ) = |«i, 2 )- Both mod¬ 
els have the same polarization vector ni j2 , which implies 
that C*i i2 = 0 for the microscopic model, too. 

In order to reveal the non-trivial topology of the mi¬ 
croscopic model, we define TR partners iTi) = (|<I> 2 ) ± 
|<f>i))/v / 2 which transform according to T|$-|-(g)) = 
|$^(— qj). Remarkably, the corresponding Chern num¬ 
bers are non-zero, C± = +2, an interesting result given 
that the general classification [18-20] rules out non-trivial 
Chern numbers for eigenstates of the Hamiltonian. The 
robustness of C± is due to a vortex like non-analyticity of 
the exciton-photon Bloch Hamiltonian Eqs. (6,7), which 
can only be removed when the splitting of the LP branch 
vanishes. As signatures of the topological phase we find 
an odd number of lines of momenta {go} encircling the 
T-point for which a 2 ({go}) = 0, experimentally acces¬ 
sible by looking for a dark excitonic state at {go}- In 
addition, we predict polaritonic edge states. Our pro¬ 
cedure of defining pairs of TR partners as suggested in 


Refs. [29, 30] is analogous to defining spin-eigenstates 
with nontrivial spin-Chern number in a QSH insulator 
with broken spin-symmetry [31, 32]. 

If TR-symmetry is broken by a small Zeeman field, 
Bz ^ 0, the pseudospin degeneracy at the TR invariant 
momenta is lifted, and two completely non-degenerate 
dispersion relations exist, see Fig. lb. For topological po¬ 
laritons with C± 7 ^ 0, we then find that for each branch of 
the dispersion the experimentally measurable pseudospin 
polarization ft(q) covers the Bloch sphere. Thus, in a TR 
broken setup, the experimentally measurable pseudospin 
polarization n(q) and the pseudospin model Eq. (1) allow 
to distinguish between topologically trivial and nontrivial 
polaritons. 

Microscopic model - The Hamiltonian of the two- 
dimensional QSH insulator in the basis of orbital states 
{| + l/2), |+3/2), |-l/2), |-3/2)} is [ 6 ] 

He(fc) = ^ e+ 0 (/c) j , H+(k) = d(k) -a, (2) 

with wavevector k and spin-orbit field d. Because of 
—* —> % 

TR-symmetry H~ (k) = H+(—k) , where a = {+,—} 
labels a pseudospin. In the following we will use the 
parametrization d x / y = Asm(k x / y ), d z = M + B(2 — 
cos (k x ) — cos(ky)) with A = Hvf/ci , B > 0, and M £ [R, 
where vf is the Fermi velocity, a the lattice spacing, and 
k is measured in a -1 . For M < 0 and B > \M\/2 the nor¬ 
malized spin-orbit field d = d/\d\ covers the Bloch sphere 
and the QSH insulator is topologically nontrivial [6] . 

An electromagnetic field A is coupled minimally to the 
semi-conductor via p —> p + eA with elementary charge 
e > 0. We work in the Coulomb gauge, linearize the 
Hamiltonian in A and expand the photon field in plane 
waves with amplitudes A qa where q is the photon mo¬ 
mentum and a labels the polarization [33]. We find for 
the optical transition matrix elements, 

k, Q) = hUd I ^ , (3) 

where /r, v label the pseudospin and band index of the 
eigenstates of Eq. (2). In order to evaluate Eq. (3) we 
use the Wigner-Eckart theorem [34] and the specific form 
of the basis states of Eq. (2) [ 6 ] . We find that optical 
transitions do not change the pseudospin a. 

A particle-hole transformation of Eq. (2) yields the 
—* —* * 

hole Hamiltonian H^(k) = —H£(—k) , with wave func¬ 
tions i(i h - = (?/) e -)*, and energies eh (k) = —e e (—k) > 0 
(the chemical potential is zero). Accounting for the at¬ 
tractive Coulomb interaction of electron and hole the 
wave function of optically active excitons is 

1*,-)=£g . (4) 

k 
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where 4>c(k) denotes the Fourier transform of the 
electron-hole wave function with respect to the relative 
coordinate. Semi-conductors typically have a large di¬ 
electric constant, which screens the Coulomb interac¬ 
tion and results in an exciton Bohr radius much larger 
as the lattice constant. Thus, the binding function 
4>c(k) is strongly peaked around k = 0 and ss 

l^?/ 2 ) ® WLq/ 2 ) with energy e x (q) « e h (q/ 2) + e e (q/ 2). 
Such excitons are described by the Hamiltonian 

H?(q) = H£(q/ 2) ® 1“ + ® H?(<?/2) , (5) 


acting on a four dimensional Hilbert space of electron- 
hole pairs. 

Diagonalizing Eq. (5) and projecting onto excitons 
gives the eigenstates \b a ) with energy e x = 2\d\. The cav¬ 
ity photon dispersion is u>q = ojq y/l + ( D/uiq ) q 2 with ujq 
set by the cavity thickness, D = h 2 Cp h /uoa 2 , photon ve¬ 
locity c p h, and momentum q measured in a . Coupling 
excitons and right (r), left (l) circularly polarized pho¬ 
tons yields polaritons. In the basis {|6 + ), |6_), | a r ), |aj)} 
the Hamiltonian takes the form 


H P 


(e x t G \ 
uil) ’ 


( 6 ) 


and the exciton-photon coupling is obtained via Eq. (3), 


go - d z )e- 2i * -(1 + 4) \ 
4 'v (1+4) -(l-4)e 2 +y 


(7) 


with 4 = d z /\d\ and e ±1<p = (4 ±id y )/\d x ± id y \. The 
coupling go oc e x (q), and is proportional to the photon 
amplitude times a numerical constant from the transi¬ 
tion matrix element Eq. (3). However, for the study of 
topological properties we can safely neglect the continu¬ 
ous (/-dependency and treat go as constant. 

Topological invariant - The Hamiltonian Eq. (6) has 
two eigenstates for both LP and UP branch: |$^ 2 ) = 
^ 12 ^ 1 , 2 ) + a i 2 l°i, 2 ) w ith 7 = {LP,UP}, exciton \bj 2 ) 
and photon pseudo-spinors \aj , 2 ). We construct new 

basis states: |$±) = (|$ 2 ) + |$7))/v^2, which are not 
eigenstates, but obey Tj$±(g)) = |<l>/j(— q)). We define 
a Chern numbers as 

C l = ~^ J , ( 8 ) 

% e BZ 

where £i :i is the Levi-Civita symbol. An explicit evalua¬ 
tion of Eq. (8) yields C^ p = +2 and C^ TP = 0 for positive 
detuning, and vice versa for negative detuning. The ex¬ 
istence of C± p = +2 is due to the vortex structure on 
the diagonal of the coupling Eq. (7). For a transition 
to a phase with different Chern numbers to occur, this 
vortex structure has to be removed, implying vanishing 
diagonal elements in Eq. (7). At the transition point the 


LP-splitting vanishes and the dispersion is twofold de¬ 
generate. Hence, the Chern numbers C± are protected 
by the splitting of LP and UP branches. 

It is enlightening to analyze the photonic component 


..lp 


) = 


,LP 


+ a 


LP 




tt r ) 


,LP 


±4 P 


e~ lv \a{) , (9) 


of |<F?jt p ). The coefficients a7 P , a 2 p G IR are continuous 
in the entire Brillouin zone, with a7 P 7 ^ 0 V q , a 2 p (r) = 
—a 7 P (r), and o 2 p (Q) = a 7 P (Q), such that a 2 p has to 
have an odd number of lines of zeros when going from the 
T-point to the Q-points. Then, the polarization vector, 
n + ~ (a pp |cr|a pp ), points north (south) at the T-point 
(Q-points) and winds twice around the 2 -axis in between. 
Thus, n± covers the Bloch sphere twice. We find that 
the polarization vector of the excitonic component |5[jl p ) 
does not cover the Bloch sphere and cannot contribute 
to Eq. ( 8 ). For all wavevectors <fo at which a 2 p = 0, 
a dark exciton eigenstate, |<I> 2 P ) = l^ 2 P )> exists. This 
constitutes a clear signature for topological polaritons. 

Breaking TR-symmetry with a Zeeman field Bz , we 
find that the LP-eigenstates |$7 P ) f° r 7 ^ 0 have 
the same topological structure as the states |4>(jl p ) for 
Bz = 0 , i.e. the |$jl p ) obtained for Bz = 0 can be con¬ 
tinuously deformed into the eigenstates |$ 7 P ) for Bz 7 ^ 0 . 
Thus, if n± covers the Bloch sphere in the case Bz = 0, 
then in the case Bz 7 ^ 0 the polarization vector 4,2 also 
covers the Bloch sphere. 

Edge states - We first evaluate the electron Hamilto¬ 
nian Eq. (2) on a cylindrical geometry with boundaries 
along the ^-direction and then couple the correspond¬ 
ing excitons to photon eigenmodes [35] . The numerically 
obtained spectrum is shown in Fig. 2a. At each bound¬ 
ary we find one pair of edge states with energies below 
the LP branch. How can these edge states be related 
to the polaritonic Chern numbers discussed above, and 
why extend the edge states all the way down to zero 
energy? Since the excitons are a direct product of topo¬ 
logically nontrivial electron and hole, it is expected to 
find non-vanishing Chern numbers = +2. The pseu¬ 
dospin degeneracy of the exciton dispersion in addition 
to TR allows to construct an operator which commutes 
with the exciton Hamiltonian and behaves exactly like 
a fermionic TR operator, for details consult our Supp. 
Mat. [35]. Due to this symmetry, the excitonic system is 
also in symmetry class All and a Z 2 index for excitons 
is given by v = \ \C+ — C* | [29]. Here, the prefactor 1/4 
was introduced because of the Chern number doubling 
due to the tensor product of hole and electron space [35]. 
We note that the coupled system of excitons and photons 
does not have any (pseudo) fermionic TR symmetry and 
is solely in symmetry class AI. The non-analyticity of the 
polariton Hamiltonian (6) does not allow to obtain a lat¬ 
tice representation and to study edge states directly. A 
lattice Hamiltonian is obtained by embedding the polari¬ 
ton space into a larger Hilbert space which contains all 
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Figure 2. (color online) Spectrum of topological polaritons for 
a system with cylindrical geometry: a) The colored dots de¬ 
pict polaritons with photonic fractions |4>j p h- The solid black 
lines show the LP and UP branches (LP and UP splitting not 
visible on this scale). As parameters A = 5|M|, B = 25\M\, 
u 0 = 3|M|/4, D = 35|M|, g 0 = |M|/24, and N = 2000 lat¬ 
tice points are used, b) The polariton branches (orange) and 
the unphysical electron-hole bands (light gray) are sketched 
schematically. Because LP-branch and negative energy bands 
carry opposite Chern numbers C± a pair of helical edge states 
(the physical part is marked purple) connects both. 

possible electron-hole (e-h) states [35], including unphys¬ 
ical e-h pairs which have either zero or negative energies, 
see Fig. 2b. Since the e-h bands at negative energy have 
Chern numbers ±2, and since the sum of all Chern num¬ 
bers in the Hilbert space of polaritons and artificial e-h 
pairs has to be zero, the polariton Chern numbers are 
compensated by the negative energy e-h Chern numbers, 
and edge states below the LP branch emerge, cf. Fig. 2. 

Effective edge model - The polaritonic edge states are 
well described by coupling excitonic edge states to pho¬ 
tonic eigenmodes [35]. To leading order perturbation the¬ 
ory we find 

e{q x ) « extez) - ’ ( 10 ) 

« HpJ ~ , (ii) 

where p = {R, L} labels right and left mover, A = ui — e x 
is the detuning between the photon Cj and exciton energy 
e x = hvp | q x |, and g the coupling strength. The excitonic 
and photonic edge-state wave functions are denoted by 
1 4’pqJ an d respectively. The right (left) moving 

exciton carries pseudospin +(—), whereas both right and 
left moving photon are linearly polarized with electric 
field parallel to the plane of incident. We find that a high 
photonic fraction of the edge state relies onj)jr~ c p h and 
go ~ wo- For e x -C 21 M our perturbative results are in 



Figure 3. (color online) Comparison of analytical (lines) and 
numerical (symbols) results for the polaritonic edge state. 
Left panel: The photon wave function (cf. Eq. (11)) is shown. 
Right panel: The coupling strength g (upper plot) and detun¬ 
ing A (lower plot) of the effective model are depicted. Same 
parameters as in Fig. 2a are used. 

very good agreement with numerics as shown in Fig. 3. 
Finally, the one-dimensional effective Hamiltonian (for 
one edge) takes the form 

He = E {e 9 b U b pq + djq a) pq a pq + (9q b l q a P q + h.C.)} . 

Q P 

( 12 ) 

The operator aJ pq creates a right, p = R, (left, p — L) 
moving photon with momentum q and dispersion oj q . In 
one dimension the elementary excitations are collective 
modes (plasmons) instead of excitons. Then, creates 
right and left moving plasmons with dispersion e q . For a 
Luttinger liquid interactions renormalize the bare Fermi 
velocity [36], and the plasmon dispersion remains linear. 

Experimental signatures - The signatures of topolog¬ 
ical polaritons are dark excitonic states along a line 
{<fo} in momentum space, and edge states below the LP 
branch, cf. Fig. la. Both are detectable via optical tech¬ 
niques [24]. Applying a Zeeman field allows to deter¬ 
mine the Chern number from analyzing the polarization 
of the bulk polaritons, see Fig. lb. Realizing topological 
polaritons in a nontrivial QSH insulator is challenging. 
However, recent developments of engineering spin-orbit 
coupling in polaritonic systems [37] and accomplishing 
slow photons in photonic crystals [38] may pave the way 
for it. 

Conclusion - We have considered a TR-invariant 
model of polaritons in a QSH cavity, and introduced a 
topological invariant for TR partners which is stabilized 
by the pseudospin splitting of polaritons. In the topolog¬ 
ical phase, polaritonic edge states below the LP branch 
exist, as well as lines in momentum space with uncoupled 
excitons, both detectable via optical techniques. 
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SUPPLEMENTAL MATERIAL 


In our main work we presented a symmetry argument that the Chern number of energy eigen¬ 
states has to be zero for the effective model. Here, we provide an argument that this result 
holds for the microscopic model, too. Furthermore, we verify that the Chern numbers of time 
reversed partners has to vanish for the effective model. Then, we show how a pseudo fermionic 
time-reversal operator for excitons can be constructed and present the embedding of the polari- 
ton Hilbert space into an enlarged Hilbert space of artificial electron-hole pair combinations. 
Finally, we consider a system with boundaries, evaluate analytically the polaritonic edge state 
wave function, and present details of the effective edge-state model. 


Chern numbers 

Microscopic model - We found that the eigenstates of 
polaritons in a QSH cavity transform under time reversal 
(TR) according to T|4>(g)) = ±|$(— q)). By symmetry 
arguments we will show that for any state which satisfies 
this behavior the corresponding Chern number defined 
via 


c = -^ J sijich* a,<i>), (s.i) 

?E BZ 

has to vanish. We decompose the integrand into 
(<0i$(«)| dj$(q)) + which is valid 

because the integration over the Brillouin zone is in¬ 
variant under inversion of momentum. Then, the 
second summand is recast: (<%<&(— q)\dj$>(—q)) = 

{d i T$(q)\d j T$(q)) = >(q)\di$(q)). We used that 

the TR-operator and partial derivative commute, and 
that T is anti-unitary. Now, Eq. (S.I) takes the form 
C ~ £ ij ({di®{$\d j ${q)) + {dj<f>{£)\di<f>(q))) = Q. 

Effective model - In our main work we have used sym¬ 
metry arguments to show that the Chern number 

C = [ d 2 q n ■ (d 9x n x d q n) . (S.2) 

47r JBZ 

of the energy eigenstate \xi, 2 ) of the TR invariant Hamil¬ 
tonian 


-Hlp = clp + 3 ■ h , (S.3) 

has to be zero. In Eq. (S.2) the polarization vector of 
the eigenstate |xi, 2 ) is = Th with h = h/\h\. We 
note that Eq. (S.2) and Eq. (S.I) are equivalent for 2x2 
Hamiltonians. For the QSH model studied in the main 
manuscript, raq 2 winds twice around the 2 -axis if q en¬ 
circles the T-point, and the effective magnetic field has 
to vanish at the T-point in order to obtain a continuous 
Hamiltonian (S.3), which results in an ill-defined polar¬ 
ization vector nq 2 - Thus, the Chern number Eq. (S.2) is 
formally not well defined. This raises the question of the 
validity of our previously presented argument. Although 
the limit q —> T is not unique, we have verified that there 
are no singular contributions from a neighbourhood of 


the T-point. Then, we can remove the T-point from the 
Brillouin zone integral, Eq. (S.2) is applicable, and our 
argument for C 12 = 0 remains valid. 

The Chern numbers C± of the TR-partners |y±) = 
(IX 2 ) ± |xi))/\/2 are defined according to Eq. (S.2) with 
polarization vector n± = (x±\®\x±)■ Now, we will 
show that C± = 0 for the effective model (S.3). To 
this end we calculate explicitly n±. For convenience we 
parametrize the effective magnetic field of Eq. (S.3) by 
h = (h±cos2(p,h.±sm2ip,h z ) T with h±(—q) = h±(q), 
ip(—q) = ip(q) + 7T and h z {—q) = —h z (q). Then, the 
energy eigenstates take the form 


_ 1 ( 'Jl T h z e 

Xl ’ 2 V2 V J 


(S.4) 


with normalized effective magnetic field h = h/\h\. Us¬ 
ing Eq. (S.4), n± can be calculated in a straightforward 
manner, 


n± = ± 


( —h z cos 2qff 
—h z sin 2 ip 

l-hi 


(S.5) 


It cannot cover the Bloch sphere, since the 2 -component 
is always positive (negative). Therefore, the Chern num¬ 
bers of the TR-partners have to vanish, C± = 0. 


Pseudo fermionic TR operator for excitons 


The exciton spectrum is degenerate with respect to the 
pseudospin a £ {+, —}, such that [H x ,a z ] = 0. Further¬ 
more, H x is invariant under TR. In a basis {|6 + ), |6_)} 
with <j z | b±) = ±|6±) the excitonic TR-operator is T = 
cr x )C, which squares to one as required by the bosonic 
statistics of excitons. The product of T and a z com¬ 
mutes with .ff x as well, and hence is a symmetry of H x , 
too. We define the anti-unitary operator Tp = <r z T. This 
is a fermionic TR-operator in the sense that Tp = ia y K. 
with Tp = — 1. 
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Enlarged polariton Hilbert space 

The BHZ Hamiltonian [6] in the basis ||+l/2), 
| + 3/2), | — 1/2), | — 3/2)} has the form 


G-(q) = — G* + (—q)a x . We emphasize that Eq. (S.9) is 
continuously defined over the entire Brillouin zone for a 
topologically nontrivial QSH insulator, and couples only 
the (physical) excitonic states to photons. 


H e ik) 


(H+(k ) 0 \ 

V 0 


(S.6) 


Polariton system on a cylindrical geometry 


where H+{k) = d{k) ■ a and H~ (k) = H+(—k) with 
pseudospin a £ {+, —}, two-dimensional wavevector k 
and spin-orbit field d. We found that optical transi¬ 
tions do not change the pseudospin a. Then, excitons 
can be characterized by a quantum number a, too. In 
order to study the topology of these excitons we em¬ 
bed the exciton Hilbert space for each pseudospin a 
into an enlarged Hilbert space; a four dimensional space 
spanned by the tensor-product states: {|o:l/2)h<8>|aT/2) e , 
|al/2) h ®|a3/2) e , |a3/2) h ® |al/2) e , |a3/2) h ® |a3/2) e }. 
Approximating H“ « Hg ® 1“ + 1“ ® yields in this 
representation 



/ 0 

d x id y 

d x + idy 

°. ^ 



* + 

II 

d x + i d y 
d x i dy 

1 

0 to 

0 

24 

d x + i dy 
d x idy 


(S.7) 


\ 0 

d x 

d x + idy 

0 ) 



Because 

of TR-symmetry 

H~{q) = 

(H+(-q))*. 

The 


Hamiltonian Eq. (S.7) describes an exciton state with en¬ 
ergy e x (q) = 2|d(<f/2)|, and three unphysical pairs with 
holes in the conduction band and/or electrons in the va¬ 
lence band. The four bands of Eq. (S.7) have Chern num¬ 
bers which result from adding the nontrivial Chern num¬ 
bers of electronic conduction and/or valence band. This 
results in a doubling of the exciton (hole-electron) Chern 
number. We find C± = =f 2 for excitons and C± = ±2 
for artificial hole-electron pairs with negative energy — e x . 
The two bands with electron and hole in the same band 
carry vanishing Chern numbers. 

The Hamiltonian of polaritons in a QSH cavity embed¬ 
ded into the extended exciton space takes the form 


H+ 

0 

G+ \ 


0 


’ 

(S.8) 

pi 

g! 

Hj 



where the exciton Hamiltonian is given in Eq. (S.7) and 
the photon Hamiltonian for right and left circularly po¬ 
larized modes is of form = ui 1 with dispersion w. 
The coupling matrix is 


G+ 


/(l —d-)dj_e 1(p -(1 + d z )d±e' v \ 
go (1 - d z ) 2 -d]_e 2iV> 

8 dle~ ^ -(1 +d z ) 2 

\(1 - d z )d±e~ lv> -(1 + d z )dj_e 1<p J 


(S.9) 


with go being a constant, d z = d z /\d\, and 

d±e 1(fi = (d x + id y )/\d\. TR-symmetry demands that 


First, we analyze the QSH insulator on a cylindrical ge¬ 
ometry with periodic boundary conditions in ^-direction 
and hard wall boundary conditions in y-direction. We 
compute the spectrum and eigenstates by numerical di- 
agonalization of the lattice Hamiltonian of Eq. (S.6). The 
eigenfunctions are plane waves in ^-direction with quan¬ 
tum numbers k x , and vanish at y = 0, L with system size 
L = L y . For given pseudospin a and wavevector k x we 
find l = 1,... ,7V conduction (valence) band eigenstates 
where L = aN with N lattice sites and lattice constant a. 
These provide a basis set for electrons \tp^ k i) an d holes 
|V’afc x z)- For a topologically nontrivial QSH insulator the 
state with l = 1 is an edge state located near a boundary. 
From now on we will label any edge state by a tilde. 

As long as momentum was a good quantum number an 
approximation for the exciton wave function was given 
by the direct product of hole and electron wave function: 

I tpaq) ~ I ^aq/ 2 ) ® W a $/ 2 ) wit h exciton momentum q. 
This approximation fails on a cylindrical geometry, since 
k y is no longer a good quantum number. Nonetheless, 
the exciton state can be expanded in product states of 
electron and hole: | i>* qxn ) = l b | wit h 

lb = leafed k x l')^Pa k x l' '•PakJ I > ( S - 10 ) 

W 


where ^ z ) are the electron and hole eigenstates in¬ 
troduced above. For sufficiently large systems we can 
approximate the excitonic wave function on a cylindri¬ 
cal geometry by projecting the excitonic eigenstates with 
periodic boundary conditions onto the electron-hole basis 
states with boundaries, namely 


I ~ ^Z c Tv{(lx) 1^/2 l C ^/2 v) . 


(S.ll) 


Cj?*(&) = <1 l>*q x /2l 




Above, the exciton state the plane wave so¬ 

lution for periodic boundary conditions with wavevec¬ 
tor q x = 2 k x , and q n = 2(2ir/L)n, where n runs from 
n = —N/2 ,..., N/2 — 1. We use the even superposition 
(plus sign) for n < 0 (cosine eigenfunction) and the odd 
one (minus sign) for n > 0 (sine eigenfunctions). 

We note that: i) The electron-hole Hilbert space is 
different for hard wall and periodic boundary conditions, 
so that lb Eq. (S.10) is the identity in the former and a 
projector in the latter space, ii) In the presence of edge 






states we replace the n = 0 state by the edge state wave with wavevector k x , spinor q (eigenvector of the a x 
function Pauli matrix with eigenvalue ±1), and real-space function 


ICJ«l^\ x /2>®fe/2> , (S.12) 

i.e. cf L ? = SiiSu'. iii) The wave function Eq. (S.ll) has 
an energy 

ex(fe,n) = E l c n"| 2 MW2,0 + e e (q x /2,l')) , (S.13) 

ll' 

where e e (eh) is the eigenvalue of the electron (hole) eigen¬ 
function of the QSH insulator on a cylindrical geometry, 
iv) We note that the expansion Eq. (S.ll) is somewhat 
analogous of projecting sine and cosine waves (free parti¬ 
cles) onto a standing wave basis (particles in a box), v) In 
the limit of infinite system size, wave function Eq. (S.ll) 
and its spectrum Eq. (S.13) converge to the solutions for 
periodic boundary conditions. 

On a cylindrical geometry the photonic eigenmodes are 


Aq x m( X >y) = \ j\ elq:DX ^ sin (9my) (S.14) 

with wavevector q x , q m = 7 r/L to , to = 1 , 2 ,..., polariza¬ 
tion vector e'er, and energy 

<*>(&:, Qm) = + , (S.15) 

with c p h as photon velocity and ujq determined by the 
thickness of the cavity. The coupling to excitons is 


iv 

X (C qx /2i'|i[^e,i] 


(S.16) 


eA 


q x m 


IC- 


g x /2// 


r lkjy) = 2^Eh— w) e \y , (S.18) 

satisfying the boundary condition i)kS.V = 0) = 0. The 
two parameters A, v (measured in inverse lattice units) 
are defined as 

v{k x ) = (S.19) 

respectively, where A,B,M are the BHZ-parameters. 
Above, all lengths (wavevectors) are measured in lat¬ 
tice units (inverse lattice units). The edge state exists 
if \k x \ < yj2\M\/B, and has an energy 

£e(k x ) = Ct hv F k x . (S.20) 

Edge states located near y = L have spinors and 
energies e^{k x ) = — ahv F k x . In the following we will 
focus on the boundary y = 0. 

Now, the coupling Eq. (S.16) is evaluated using the 
approximation Eq. (S.12) and the result Eq. (S.17). We 
find that only p-polarized light (linearly polarized light 
with electric field parallel to the plane of incident which is 
perpendicular to the y-direction) couples: g“ CT=p / 

0 , whereas s-polarized light (electric field perpendicular 
to the boundary) does not: g^ =s = 0. We find 

5m(M = IVK(y) 2 M™y )> (S-2i) 

for both pseudospins a = ±. 

As long as the exciton bulk states are energetically 
much higher as the edge state, e x = fov F \q x \ <C 2|M|, we 
can neglect those and use an effective description, 


where c, v label the conduction and valence band, respec¬ 
tively. 

The polaritonic modes are obtained by evaluating the 
coupling Eq. (S.16) numerically and diagonalizing the 
corresponding polariton Hamiltonian with exciton en¬ 
ergy Eq. (S.13) and photon energy Eq. (S.15) for given 
wavevector q x and pseudospin a. We find pairs of po¬ 
laritonic edge states lying energetically below the lower 
polariton branch. 


Polaritonic edge state model 

Replacing k y —> —id y in Eq. (S.6) allows us to calculate 
analytically the edge-state wave function for hard wall 
boundary conditions in y-direction, see Ref. [2, 39, 40]. 
For wave functions located near y = 0 [41] we find 

Csjg) = ? 7fc x (y) </>+i > (S- 17 ) 


/ exfe) 0 g + (qx)\ 

H' F (q x )= 0 Uqx) g-(qx) , (S. 22 ) 

\(5 + (fe)) t (5"(w)) t u(q x ) J 

where g a is a row vector in the photon-space, see 
Eq. (S.21), and u(q x ) = u}(q x ,q m )S mm f a diagonal ma¬ 
trix, see Eq. (S.15). In leading order degenerate pertur¬ 
bation theory we find 


e(&) ~ £(o ){qx) + £( 2 )(5*) 

r \ _ |gm(ga/2)| 2 

M ^ A m (q x ) 

«1%) +1*&> 

( g ±(g a /2))* 
( Qx ) 


= 1 


(S.23) 


\AZ P ) . (S.24) 


Above, we have defined the detuning A m (q x ) = 
w(q x , qm) — £x(g£c), and an index p which labels right (R) 
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and left (L) moving edge states. The excitonic compo¬ 
nent of the right (left) mover has pseudospin a = +(—), 
whereas the photonic component is p-polarized indepen¬ 
dently of the direction of motion. Since is equal for 
both pseudospins, the edge-state energy is degenerate for 
p = {R,L}. Evaluating the photonic part of & pq:c (y) 
Eq. (S.24) results in a wave function which is exponen¬ 
tially localized near the boundary, too. We define the 

photon fraction as F = \J |P p h | 2 where P p h projects 

onto the photonic part, so that Eq. (S.24) yields 


ergy of the localized photon wave function and g an ef¬ 
fective coupling. Diagonalizing Eq. (S.26) and expanding 
in powers of g /A yields 

KQx) ~ e x (q x ) - (f Afe) , (S.27) 

\* Pq J « \r P j - > (s.28) 

^\Qx) 

with detuning A = Q — e x . This allows us to extract 


F(q x ) * 


\g™{q x /2)\' 2 

A m (g x ) 2 


(S.25) 


A simple effective model for right or left moving po- 
laritonic edge states is 

=((MSr StO ■ (S - 26) 

where e x (q x ) = hvp\q x is the exciton energy, Co the en- 




1 v- l^fe/ 2 )] 2 


F{q x ) m A m {qx) 

Xf ^ lgm(ga :/ 2 )| 2 

[qx) F{q x f^ A m (q x ) ’ 

, Tph v = 1 (ffm(W 2 ))* I Z*=P\ 

\F P q J F{qx) 2 L A m (q x ) 1 ’ 

by comparing with Eq. (S.23) and Eq. (S.24). 


(5.29) 

(5.30) 

(5.31) 













